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^ ■ An updated version is presented of the Fortran code FeynHiggs for the evaluation 

of the neutral CP-even Higgs sector masses and mixing angles. It differs from the 
previous version by a modification of the renormalization scheme concerning the treat- 
ment of subleading terms at the one-loop level; the two-loop corrections, for which the 
leading contributions of (D(ata s ) and C(a 2 ) are implemented, are not affected by the 
modified renormalization prescription. Besides stabilizing the Higgs mass calculations 
and avoiding unphysically large threshold effects, the mass of the lightest MSSM Higgs 
boson, m/j, is increased by 1-2 GeV for most parts of the MSSM parameter space. 
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Abstract 

An updated version is presented of the Fortran code FeynHiggs for the evaluation 
of the neutral CP-even Higgs sector masses and mixing angles. It differs from the 
previous version by a modification of the renormalization scheme concerning the treat- 
ment of subleading terms at the one-loop level; the two-loop corrections, for which the 
leading contributions of 0(ata s ) and 0(a^) are implemented, are not affected by the 
modified renormalization prescription. Besides stabilizing the Higgs mass calculations 
and avoiding unphysically large threshold effects, the mass of the lightest MSSM Higgs 
boson, mh, is increased by 1-2 GeV for most parts of the MSSM parameter space. 

1 Introduction 

The search for the lightest Higgs boson in the Minimal Supersymmetric Standard Model 
(MSSM) is one of the main goals at the present and the next generation of colliders. Therefore 
the precise knowledge of the dependence of the masses and mixing angles of the Higgs sector 
of the MSSM on the relevant supersymmetric parameters is of high importance. 

In this note we present an updated version of the Fortran code FeynHiggs that eval- 
uates the neutral CP-even Higgs sector masses and mixing angles || [J. It differs from the 
previous version as presented in Ref. by a modification of the renormalization scheme 
concerning the treatment of subleading terms at the one-loop level; the two-loop corrections, 
for which the leading contributions of 0(a t a s ) and 0(of ) are implemented, are not affected. 
In particular, an MS renormalization for tan f3 and the field renormalization constants has 
been used (where the MS quantities are evaluated at the scale m t ). The renormalization in 
the new version of FeynHiggs does no longer involve the derivative of the A boson self-energy 
and the AZ mixing self-energy. This leads to a more stable behavior around thresholds, e.g. 
at Ma ~ 2m t , and avoids unphysically large contributions in certain regions of the MSSM 
parameter space. Thus, the new renormalization scheme stabilizes the prediction of the 
masses and mixing angles in the CP-even Higgs sector of the MSSM. 

2 Renormalization schemes 

The Higgs sector of the MSSM |4| consists of two neutral CP-even Higgs bosons, h and 
H (m h < m H ), the CP-odd A boson, and two charged Higgs bosons, H ± . At the tree- 
level, m h and m H can be evaluated in terms of the Standard Model (SM) gauge couplings 
and two additional MSSM parameters, conventionally chosen as Ma and tan ft, the ratio of 
the two vacuum expectation values (tan/5 = i^/fi)- Beyond lowest order, the Feynman- 
diagrammatic (FD) approach allows to obtain in principle the most precise evaluation of the 
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neutral CP-even Higgs boson sector, since in this way the effect of different mass scales of 
the supersymmetric particles and of the external momentum can consistently be included. 
The masses of the two CP-even Higgs bosons are obtained in this approach by determining 
the poles of the h — //-propagator matrix, which is equivalent to solving the equation 

[ q 2 - Ktree + Mf) } [ f ~ ^l.tree + ] " [ t hH {q 2 ) f = 0, (1) 

where £ s ,s = h,H,hH, denote the renormalized Higgs boson self-energies. For the renor- 
malization within the FD approach usually the on-shell scheme is applied ||. This means in 
particular that all the masses in the FD result are the physical ones, i.e. they correspond to 
physical observables. Since eq. (|I|) is solved iteratively, the result for m h and m# contains a 
dependence on the field renormalization constants of h and H, which is formally of higher 
order. Accordingly, there is some freedom in choosing appropriate renormalization condi- 
tions for fixing the field renormalization constants (this can also be interpreted as affecting 
the renormalization of tan/3). Different renormalization conditions have been considered, 
e.g. (£' denotes the derivative with respect to the squared momentum): 

1. on-shell renormalization for T, z , ^U, S^, Y<az, and 5v\/vi = 5v 2 /v2 [§] 

2. on-shell renormalization for T, z , E^, XUz, and 5vi = 5v it di v , i — 1,2 

3. on-shell renormalization for @, MS renormalization for 5Z^,5Zh, tan/5 [[?[]. 

The previous version of FeynHiggs is based on renormalization 1 , involving the derivative of 
the A boson self-energy. The new version of FeynHiggs, see www.feynhiggs.de, is based on 
renormalization 3 (a detailed discussion can be found in Ref. [[?]]). 

3 Numerical comparison 

In this section we numerically compare the output of the previous version (based on renor- 
malization 1) and the new version (based on renormalization 3) of FeynHiggs. We also show 
results for the recently obtained non- logarithmic 0{aj) corrections |5], [J] that are also in- 
cluded in the new version of FeynHiggs. The comparison is performed for the parameters of 
the three LEP benchmark scenarios [|10J. In this way, the effect of the new renormalization 
and the non-logarithmic O(of) corrections on the analysis of the LEP Higgs-boson searches 
can easily be read off. 

In Figs. [I|-|3] we show the results in the "m™ ax " , "no-mixing" and "large /i" scenario 
as a function of Ma (left column) and of tan/3 (right column) for two values of tan/3 
(tan (3 = 3, 50) and Ma {Ma = 100, 1000 GeV for the m™ ax and the no-mixing scenario, 
Ma = 100,400 GeV for the large /i scenario), respectively. The solid lines correspond to the 
new result while the dashed lines show the old results. The dotted lines correspond to the 
new result including the non-logarithmic 0(af) contributions. Concerning the new renor- 
malization scheme, in the m™ ax (Fig. HD and the no-mixing scenario (Fig. |2|) the new result 
is larger by GeV for not too small Ma and tan/3. For small tan/3 and large Ma the en- 
hancement can be even larger. In the large /i scenario (Fig. the largest deviations appear 
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for small tan/3 for both large and small Ma- While the previous prescription for the field 
renormalization constants leads to unphysically large threshold effects in some regions of the 
parameter space, which arise from the AZ mixing self-energy and the derivative of the A bo- 
son self-energy, no threshold kinks are visible for the result based on the new renormalization. 
The shift in of ~l-2 GeV related to the modification of the renormalization prescription 
lies in the range of the anticipated theoretical uncertainty from unknown non-leading elec- 
troweak two-loop corrections W% - The new 0(af) corrections can further increase m/, by up 
to ~3 GeV for large t mixing (a detailed analysis will be presented elsewhere |p~2f ). 

The new version of FeynHiggs can be obtained from www.feynhiggs.de . 
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Figure 1: The new renormalization (3, solid) and the old scheme (1, dashed) are compared 
in the m™ ax scenario. The dotted line shows the inclusion of the non-logarithmic O(a^) 
corrections. The lower curves are for tan/5 = 3 (left plot) or M A = 100 GeV (right). The 
upper curves are for tan/3 = 50 (left) or M A = 1000 GeV (right). 




Figure 2: The new renormalization (3, solid) and the old scheme (1, dashed) are compared 
in the no-mixing scenario. The dotted line shows the inclusion of the non-logarithmic O(af) 
corrections. The lower curves are for tan/3 = 3 (left plot) or M A = 100 GeV (right). The 
upper curves are for tan/3 = 50 (left) or M A = 1000 GeV (right). 
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Figure 3: The new renormalization (3, solid) and the old scheme (1, dashed) are compared 
in the large \x scenario. The dotted line shows the inclusion of the non- logarithmic 0(af) 
corrections. The lower curves are for tan/3 = 50 (left plot) or M A = 100 GeV (right). The 
upper curves are for tan/3 = 3 (left) or M A = 400 GeV (right). 
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